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In  this  article  we  address  two  questions:  Why  do  freely  evolving  vortices  weaken  on  average, 
even  when  the  viscosity  is  very  small?  Why,  in  the  fluid’s  interior,  away  from  vertical  boundaries 
and  under  the  influence  of  Earth’s  rotation  and  stable  density  stratification,  do  anticyclonic 
vortices  become  dominant  over  cyclonic  ones  when  the  Rossby  number  and  deformation 
radius  are  finite?  The  context  for  answering  these  questions  is  a  rotating,  conservative, 
Shallow-water  model  with  Asymmetric  and  Gradient-wind  Balance  approximations.  The 
controlling  mechanisms  are  vortex  weakening  under  straining  deformation  (with  a  weakening 
that  is  substantially  greater  for  strong  cyclones  than  strong  anticyclones)  followed  by  a  partially 
compensating  vortex  strengthening  during  a  relaxation  phase  dominated  by  Vortex 
Rossby  Waves  (VRWs)  and  their  eddy-mean  interaction  with  the  vortex.  The  outcome  is  a 
net,  strain-induced  vortex  weakening  that  is  greater  for  cyclones  than  anticyclones  when  the 
deformation  radius  is  not  large  compared  to  the  vortex  radius  and  the  Rossby  number  is 
not  small.  Furthermore,  when  the  exterior  strain  flow  is  sustained,  the  vortex  changes  also 
are  sustained:  for  small  Rossby  number  (i.e.,  the  quasigeostrophic  limit,  QG),  vortices  continue 
to  weaken  at  a  relatively  modest  rate,  but  for  larger  Rossby  number,  cyclones  weaken  strongly 
and  anticyclones  actually  strengthen  systematically  when  the  deformation  radius  is  comparable 
to  the  vortex  radius.  The  sustained  vortex  changes  are  associated  with  strain-induced  VRWs  on 
the  periphery  of  the  mean  vortex.  It  therefore  seems  likely  that,  in  a  complex  flow  with  many 
vortices,  anticyclonic  dominance  develops  over  a  sequence  of  transient  mutual  straining  events 
due  to  the  greater  robustness  of  anticyclones  (and  occasionally  their  net  strengthening). 
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1.  Introduction 

Observational  and  computational  studies  have  shown  that  coherent  vortices  emerge 
in  large-scale  geophysical  fluid  flows  under  many  conditions,  in  particular  —  with 
relevance  to  the  present  study  —  in  the  fluid’s  interior  away  from  solid  boundaries 
and  under  the  influences  of  planetary  rotation  and  stable  density  stratification. 
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Terrestrial  examples  include  atmospheric  vortices  near  the  tropopause  and  oceanic 
vortices  in  the  thermocline  and  abyss.  Other  examples  can  be  found  on  the  giant 
planets,  most  spectacularly  the  Great  Red  Spot  of  Jupiter  (figure  1).  Coherent 
vortices  have  also  been  found  in  computational  turbulence  studies  of  initially 
randomly  distributed  vorticity  in  stratified,  rotating  environments  (McWilliams 
et  al.  1994)  where  axisymmetrization,  vertical  alignment,  and  successive  mergers  of 
like-sign  vortices  cause  the  aggregation  of  vorticity  into  a  decreasing  number  of 
larger  vortices. 

The  emergent  vortices  tend  to  be  relatively  axisymmetric  around  a  core  vorticity 
extremum.  This  is  true  in  part  because  an  axisymmetric  azimuthal  circulation  is  a 
stationary  solution  to  the  conservative  fluid  equations  in  gradient-wind  and  hydrostatic 
balance,  and  further  because  linearly  stable  vortices  are  often  robust  to  even  substantial 
perturbations  that  deform  them  away  from  their  axisymmetric  stationary  state.  When 
perturbed,  vortices  tend  to  relax  back  to  an  axisymmetric  state  on  an  advective  time 
scale,  exhibiting  a  long-time  weakening  after  many  disturbance  and  relaxation  cycles. 
The  final  axisymmetrization  is  a  result  of  both  advection  and  diffusion,  but  diffusion 
alone  is  insufficient  to  account  for  the  degree  of  weakening.  This  relaxation  process 
is  often  referred  to  as  vortex  axisymmetrization,  and  many  studies  have  tried  to  capture 
its  underlying  dynamics.  Melander  et  al.  (1987)  developed  a  descriptive  framework 
to  explain  the  vortex  axisymmetrization  process.  They  showed  that  elliptical  vortices 
have  nonzero  potential  vorticity  located  outside  the  saddle  points  in  the  co-rotating 
streamfunction,  causing  the  potential  vorticity  to  be  pulled  away  from  the  vortex  and 
filamented,  leaving  behind  a  smaller  circular  vortex.  Montgomery  and  Kallenbach 
(1997)  showed  that  vorticity  perturbations  on  smooth-shaped,  monopolar  vortices 
propagate  outward  as  Vortex  Rossby  Waves  (VRWs)  that  are  azimuthally  sheared 


Figure  1.  Jupiter’s  long-lived  anticyclonic  vortices  visible  as  pale  ovals.  This  photograph  is  from  Voyager  I 
on  13  February  1979.  The  moons  Io  (left)  and  Europa  are  also  visible,  and  the  Great  Red  Spot  vortex  lies 
behind  Io  (NASA/JPL). 
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by  the  differential  angular  velocity  of  the  mean  axisymmetric  vortex.  These  fluctuations 
eventually  transfer  their  energy  back  to  the  mean  vortex,  thereby  strengthening 
it  relative  to  the  initial  azimuthal  mean  flow.  In  a  complementary  study,  Schecter 
et  al.  (2000)  demonstrated  some  circumstances  where  energy  exchange  occurs  between 
a  fluctuation,  sometimes  referred  to  as  a  quasi-mode  (Briggs  et  al.  1970),  and  the  mean 
vortex  at  a  critical  radius  where  the  azimuthal  wave  speed  of  the  fluctuation  equals 
the  angular  speed  of  the  vortex. 

To  further  understand  the  robustness  of  vortices,  Bassom  and  Gilbert  (1999; 
hereafter  BG99),  examined  how  a  barotropic  vortex,  with  a  smooth  radial  vorticity 
profile  responds  to  an  external  strain  flow  that  is  impulsive  in  time,  has  an  extended 
duration,  or  is  somewhat  random.  They  analyzed  the  weakening  of  the  azimuthal- 
mean  vortex  during  the  straining  period  and  the  subsequent  relaxation  of  the  vortex 
back  to  an  axisymmetric  state.  They  found  that  the  net  effect  of  an  external  strain 
flow  is  a  weakening  of  the  original  vortex. 

In  this  article  we  further  explore  this  phenomenon  using  the  Shallow-Water 
Equations  (SWEs)  with  small  but  finite  Rossby  and  Froude  numbers  (i.e., 
Ro  =  V/fL,  where  V  and  L  are  vortex  azimuthal  velocity  and  radial  length  scales, 
and /is  the  Coriolis  frequency;  Fr  =  V/^fgH,  where  g  is  the  gravitational  acceleration, 
and  H  is  the  resting  depth  of  the  fluid  layer).  A  striking  characteristic  of  coherent  vor¬ 
tices  —  not  yet  well  understood  —  is  the  fact  that,  with  finite  Ro  and  deformation 
radius  (i.e.,  Ld  —  gH/f ),  anticyclones  tend  to  be  more  prevalent  than  cyclones  in 
flows  that  are  not  strongly  influenced  by  horizontal  or  vertical  boundaries.  (This  is  in 
contrast  with  the  development  of  stronger  cyclones  in  baroclinic  flows  dominated  by 
boundary  potential  temperature  gradients,  e.g.,  synoptic  winter  cyclones  in  the  atmo¬ 
sphere.)  Many  studies  have  noted  the  preferential  existence  of  anticyclones  both  in 
nature  (e.g.,  atmospheric  blocking  anticylones  and  submesoscale  coherent  vortices  in 
the  oceanic  interior;  McWilliams  1985)  and  in  computational  solutions  of  the  SWE 
or  its  balanced  approximation  (Cushman-Roisin  and  Tang  1990,  Arai  and  Yamagata 
1994,  Polvani  et  al.  1994,  Cho  and  Polvani  1996)  or  in  a  fully  three-dimensional  (3D) 
flow  (Yavneh  et  al.  1997).  Whether  there  is  a  discrepancy  in  the  number  of 
emerging  vortices  with  anticyclonic  vorticity  versus  cyclonic  vorticity,  whether  cyclones 
are  more  susceptible  to  an  instability,  or  whether  anticyclones  are  more  robust  to 
externally  induced  disturbances  are  all  issues  of  continuing  investigation. 

Cushman-Roisin  and  Tang  (1990)  showed  a  preferential  emergence  of  anticyclones 
from  a  random  vorticity  field  when  characteristic  sizes  were  larger  than  Ld.  They 
suggested  this  was  due  to  an  inherent  tendency  of  cyclones  to  self-destruct.  They 
proposed  that  the  quadratic  component  of  the  planetary  /1-effect  in  their  model  is 
responsible  for  balancing  Rossby-wave  dispersion  in  anticyclones,  akin  to  a  soliton 
solution  (cf  Matsura  and  Yamagata  1982).  Cyclones,  on  the  other  hand,  would  fail 
to  remain  coherent  since  no  such  dynamical  balance  is  possible.  They  also  advanced 
a  heuristic  stability  argument  favoring  anticyclone  survival.  However,  since  anticyclonic 
dominance  has  also  been  observed  on  an  /plane  (e.g.,  Polvani  et  al.  1994),  this 
/1-dispersion  argument  cannot  be  the  generic  explanation. 

Arai  and  Yamagata  (1994)  also  observed  anticyclonic  dominance  in  numerically 
simulated  flows  when  Ro  is  sufficiently  large.  To  explain  the  phenomenon,  they 
examined  the  problem  of  elliptical  vortices  on  an  /-plane.  Anticyclonic  ellipses 
were  observed  to  undergo  the  familiar  axisymmetrization  process,  regardless  of  the 
magnitude  of  Ro  or  LdjL.  In  contrast,  cyclonic  ellipses  elongated  and  eventually  split 
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into  two  vortices.  This  tendency  to  elongate  and  split  decreases  as  the  ratio  LdfL 
decreases  (see  their  figure  8).  They  suggested  that  this  is  the  mechanism  for  anticyclonic 
dominance  and  remarked  that  vortex  splitting  is  associated  with  the  development 
of  additional  saddle  points  in  the  co-rotating  streamfunction  near  the  vortex  core 
that  implies  an  advective  tendency  to  pull  the  vortex  apart.  We  show  that  anticyclonic 
dominance  does  not  require  vortex  splitting.  To  identify  flow  regimes  displaying 
anticyclonic  dominance,  Polvani  et  al.  (1994)  performed  numerical  simulations  of  the 
SWE  for  the  unforced  evolution  of  initially  random  vorticity  fields  with  a  wide  range 
of  Ro  and  Fr  values.  They  found  that  as  Fr  (—  Ro  L/Ld)  increases  at  fixed  Ro, 
so  does  the  negative  value  of  vorticity  skewness,  implying  increased  anticyclonic 
dominance.  This  is  in  apparent  contradiction  with  the  tendency  favoring  the  occurrence 
of  vortex  splitting  (Arai  and  Yamagata  1994).  In  order  to  reconcile  this  discrepancy, 
either  initial  condition  differences  are  important  (e.g.,  whether  or  not  the  vortices 
have  zero  circulation)  or  the  interpretation  that  cyclonic  elongation  and  splitting  are 
the  unique  cause  of  anticyclonic  dominance  is  incorrect. 

Finear  stability  studies  for  vortices  have,  for  the  most  part,  failed  to  find  a  greater 
instability  for  cyclones.  There  are  two  finite-Ro  instabilities,  centrifugal  (Charney 
1973,  Floskins  1974)  and  ageostrophic  anticyclonic  (Molemaker  et  al.  2004),  that 
imply  greater  anticyclonic  instability.  An  exception  is  the  study  by  Stegner  and 
Dritschel  (2000)  that  considered  the  instability  and  finite  amplitude-evolution  of 
zero-circulation  vortices.  Although  they  show  a  greater  cyclonic  growth  rate  of 
instability  compared  to  anticyclones,  the  instability  is  essentially  a  quasi-geostrophic, 
inflection-point  type,  with  only  modest  differences  in  evolution  arising  between 
cyclones  and  anticyclones  with  increasing  Ro  and  increasing  sharpness  of  the  radial 
vorticity  profiles.  Therefore,  we  believe  that  linear  stability  behaviors  are  an  insufficient 
explanation  for  the  preferential  persistence  of  anticyclones  observed  in  interior  flow 
regimes  where  many,  if  not  most,  vortices  have  nonzero  circulation.  Many  different 
radial  profile  shapes  arise  for  vortices  in  complex  flows,  and  by  their  survival  they 
are  evidently  stable. 

Cho  and  Polvani  (1996)  showed  that  geostrophic  adjustment  by  inertia-gravity  wave 
radiation  in  the  SWE  favors  stronger  end-state  anticyclones  at  intermediate  Ro  values, 
and  Yavneh  and  McWilliams  (1994)  showed  that,  in  the  core  of  strongly  deformed 
cyclones  when  the  fluid  depth  approaches  zero  in  the  SWE,  the  dynamical  approxima¬ 
tion  of  diagnostic  force  balance  breaks  down,  which  instigates  an  enhanced  dissipation 
rate.  While  both  of  these  processes  may  contribute  to  anticyclonic  dominance  in 
general,  neither  is  relevant  to  the  essential  phenomenon  that  occurs  in  balanced 
vortex  dynamics. 

In  this  article  we  present  a  new  theory  for  the  cause  of  anticyclonic  dominance  in 
large-scale  geophysical  flows,  based  on  the  demonstration  that  anticyclones  are 
more  robust  than  cyclones  to  perturbative  disturbances  at  finite  Ro.  The  theory’s 
ingredients  are  the  diagnostic  force  balances  for  the  azimuthally  averaged  vortex 
and  its  asymmetric  fluctuations  ( balanced  dynamics)',  a  linearized  approximation  for 
the  fluctuation  dynamics  in  the  presence  of  the  mean  vortex;  and  a  fully  nonlinear 
representation  of  the  azimuthally  averaged  momentum  and  mass  fluxes  by  the 
fluctuations  (quasi-line ar,  eddy-mean  interaction',  e.g.,  Findzen  and  Holton  1968). 
Thus,  the  theory  does  not  consider  extreme  finite-amplitude  events,  such  as  vortex 
splitting  and  filament  expulsion  (e.g.,  Nof  1990,  1991,  Arai  and  Yamagata  1994, 
Drijfhout  2003)  that  we  find  not  to  be  essential  for  the  generic  evolutionary  behavior. 
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One  way  to  describe  the  disturbance  of  any  particular  vortex  in  a  complex,  turbulent 
flow  is  to  view  the  external  influences  as  a  time-varying  strain  flow  acting  on 
that  vortexf.  As  shown  in  BG99,  a  strained  barotropic  vortex  undergoes  an  initial 
weakening  followed  by  an  incomplete  recovery  to  axisymmetry,  resulting  in  a  net 
weakening.  We  show  that  the  essential  difference  between  cyclones  and  anticyclones 
is  in  their  response  to  such  an  external  strain,  whereby  cyclones  exhibit  greater  net 
weakening  J. 

Section  2  is  the  model  derivation  for  the  influence  of  an  external  strain  flow  on  a 
vortex  initially  in  axisymmetric,  gradient-wind,  and  hydrostatic  balance.  Section  3 
poses  our  specific  problems  and  describes  the  numerical  methods  employed. 
Section  4  is  a  preview  of  the  essential  phenomena.  Section  5  analyzes  the  weakening 
of  a  vortex  during  the  impulsive  straining  period.  Section  6  describes  the  strengthening 
of  a  vortex  as  it  recovers  from  the  external  strain  through  VRW  evolution.  Section  7 
presents  the  sensitivities  to  different  parameters,  such  as  strain  duration,  Ro  and  Fr, 
and  advective  nonlinearity  of  the  fluctuations.  Section  8  examines  solutions  when  the 
external  strain  flow  is  sustained  in  time.  Finally,  section  9  provides  a  summary  of 
the  study  and  a  discussion  of  remaining  issues. 


2.  Dynamical  model 

We  investigate  the  effect  of  a  specified  external,  large-scale,  straining  flow  on  an  initially 
axisymmetric,  balanced  vortex.  The  strain  excites  an  asymmetric  perturbation  in  the 
vortex  that  subsequently  relaxes  towards  axisymmetry  after  the  strain  abates.  We  use 
the  conservative,  rotating  SWE  with  a  free  upper  surface  and  flat,  solid  bottom  surface. 
In  cylindrical  coordinates  (r,  X,  t),  the  primary  dependent  variables  are  the  radial 
and  azimuthal  velocity  ( u ,  v)  and  geopotential  </>.  In  a  companion  article  (McWilliams 
et  al.  2003;  hereafter  MGM03),  the  model  equations  are  fully  developed  in  the 
context  of  VRW  solutions  and  their  interaction  with  the  mean  vortex.  Here,  more 
briefly  we  record  the  relevant  equations  extended  to  include  an  external  strain  flow. 

We  decompose  the  primary  variables  with  an  ordering  parameter,  e,  that  measures 
the  relative  sizes  of  the  vortex  fluctuation  and  external  strain  rate  compared  to  the 
undisturbed  vortex: 


u  —  cm-  +  e2(Ui)  +  eu'e  +  e2{ue) 

v  =  v(r)  +  evl  +  e2  <vf)  +  ev'e  +  e2  (ve}  (1) 

(p  =  cp(r )  +  etp'j  +  e2  {(pi)  +  €(p'e  +  e2{cpe). 


fEarlier  work  by  Brickman  and  Ruddick  (1990)  examined  the  instability  of  anticyclonic  vortices  to  a 
barotropic  straining  flow.  While  the  setup  is  similar  to  the  problem  solved  here,  they  did  not  consider  the 
comparative  cyclonic  case  as  we  do  here,  and  their  model  was  restricted  to  vortices  with  a  discontinuous 
potential-vorticity  distribution  (i.e.,  an  anticyclonic  ‘lens’)  rather  than  the  continuous  profiles  analyzed  here. 
JFully  nonlinear  simulations  in  3D  by  Graves  (2005)  examining  free-decay  turbulence,  vortex  axisymmetriza- 
tion,  and  vortex  merger  at  finite  Rossby  number  strongly  support  the  paradigm  of  strain-weakening/ 
relaxation-strengthening/anticyclonic-dominance  developed  here  by  our  semi-analytical  analysis. 
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The  overbar  denotes  the  initial  vortex  that  is  a  stationary  state  in  the  absence  of  external 
strain.  The  primes  and  brackets  denote  azimuthally  asymmetric  and  symmetric  parts  of 
the  strain-induced  flow.  The  subscripts  e  and  i  refer  to  the  ‘external’  strain  flow  and 
‘internal’  vortex  perturbation,  respectively.  The  and  •'  quantities  are  referred  to 
as  the  vortex  changes  and  fluctuations.  Note  that  we  have  subtracted  off  any  bulk 
external  velocity  (e.g.,  evaluated  at  the  vortex  center,  r  =  0),  which  only  induces  a 
bulk  displacement  of  the  vortex,  to  focus  on  the  vortex  deformation  induced  by  the 
external  velocity  gradients  across  the  vortex.  The  azimuthally  averaged  part  of 
the  external  strain  flow  is  an  order  smaller  than  the  asymmetric  part. 

Inserting  (1)  into  the  SWE  and  separating  by  powers  in  e,  we  obtain  at  O(e0)  the 
stationary  gradient-wind  balance  for  v  and  </>, 


d(f> 

dr 


=fy  + 


v2 


(2) 


and  at  0(e)  a  linear  SWE  system  for  the  fluctuations, 

Du'i  -  fv'  =  —3, .<//.  +  n(2Ve  -  3;X) 

Dv\  +  r)ui  =  -  i  3;0'  -  Cm',  -  f!3x v'e  (3) 

Drp'j  +  ghS'j  +  u'j  d,4>  —  —ghS'e  -tod  X(p'e  -  u'e  dr<j>. 


There  are  no  terms  related  to  the  external  flow  by  itself  since  it  is  assumed  to  satisfy  its 
own  dynamical  balance  equations  in  the  absence  of  the  vortex.  Relevant  quantities  are 
defined  by 


D-d,  +-dx  ,  Q  —  ,  C  =  9,0‘v), 
r  r  r 

b—  /+?>  f=/+2£2,  h  —  H+~,  (4) 

g  v 

-  9  -2  hi  1  a  .  1  a  , 

q  =  =,  Y  9  =-dr(ru)  +  -dxV, 

h  gh  r  r 

where  Q  is  the  angular  velocity  of  the  vortex;  D  is  the  substantial  time  derivative  asso¬ 
ciated  with  advection  by  the  mean  vortex;  £  is  its  relative  vorticity;  rj  is  its  absolute 
vorticity;  f  is  twice  its  absolute  rotation  rate;  h  is  its  layer  depth;  q  is  its  potential 
vorticity;  y~ 1  is  its  local  deformation  radius  (~  Lj)\  and  S'  is  the  asymmetric  component 
of  horizontal  divergence. 

To  exclude  the  fast  inertia-gravity  wave  oscillations  in  (3),  we  make  an  approxi¬ 
mation  based  on  diagnostic  force  balance  (McWilliams  2003),  specifically  the 
Asymmetric  Balance  (AB)  approximation  of  Shapiro  and  Montgomery  (1993) 
appropriate  to  nearly  axisymmetric  baroclinic  flowsf  (assured  here  by  e«l).  AB 
involves  an  iterative  substitution  for  the  velocity  in  the  momentum  equations  in  (3), 


tThe  accuracy  of  AB  solutions  relative  to  the  hydrostatic  Primitive  Equations  for  vortex  flows  has  been  well 
tested  (e.g.,  by  Moller  and  Montgomery  (1999)  and  Schecter  and  Montgomery  (2003)  for  vortices  with 
finite  Ld).  AB  becomes  inaccurate  in  the  barotropic  limit  at  finite  Ro,  and  we  avoid  its  use  in  this  situation. 
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starting  from  a  generalized  geostrophic  balance  in  the  first  iterate  (i.e.,  equating 
the  second  left-side  and  first  right-side  terms  in  the  first  two  equations).  After  the 
next  iteration,  the  result  is 


M,  =  -  - 


“  9x0;  +  =-£>(9r0;)  +  ^U.'e  +  f29yV^, 


V;  H 


il  _r 

9,-0'  -  -£(9x0;-)  +  £2(9x4  -  20 

rjr 


(5) 


Substituting  (5)  into  the  continuity  equation!  in  (3)  and  multiplying  everything  by  — y2, 
we  get  the  fluctuation  geopotential  tendency  equation  (cf  MGM03): 


D[qah]~  X-  Gdkti=-m 


+  r 


1  „  (r 


+  ^4))  +  ^ (29xv'e  -  dlu'e)  -  ge 
—ifu’e  -  mye)  + 

?y 


(6) 


where 


A  =  V20'  -  y20',  G(r)  =  £9,M,  V2  =  ^9,(^9,)  +  4  9^-  (7) 

r  \y-  /  r- 

is  the  fluctuation  potential  vorticity  in  the  AB  approximation.  All  fluctuation 
(primes)  and  vortex  change  terms  (brackets)  are  initially  zero,  as  is  the  external  strain 
(section  3).  The  boundary  conditions  for  0-,  and  thus  for  9,0'.,  are  zero  for  both 
itself  and  its  radial  derivative  at  r  =  0,  and  decay  ~r-2  at  the  outer  boundary  due  to 
the  choice  for  the  external  strain  flow  (section  3)  or  even  more  steeply  ~  e~rlLd  after 
the  strain  abates. 

Inserting  (1)  into  the  SWE  and  azimuthally  averaging,  we  obtain  at  0(e2)  beyond 
(2)  the  following  system  for  the  vortex  changes: 

9r(«;>  -  l<V;>  =  9r(0;)  -  (EAO  +  2H(ve) 

9;(V;>  +  h(M;>  =  —  (£ArV)  -  l{ue)  (8) 

9f(0i>  +-dr(rgh(Ui))  -  -(XN*)  —  - dr(rgh(ue)), 

r  r 


fAs  discussed  in  MGM03,  there  is  an  alternative  derivation  for  the  AB  asymmetric  geopotential  tendency 
equation,  where  the  velocity  approximations  (5)  are  substituted  into  the  linearized  asymmetric  potential- 
vorticity  equation.  The  two  alternatives  for  the  AB  approximation  are  referred  to  as  the  continuity  path 
and  potential-vorticity  path,  respectively.  The  paths  are  of  an  equivalent  formal  order  of  accuracy  in  Ro,  with 
errors  0((\D\/P)2 ),  but  they  differ  in  several  of  their  higher-order  terms.  We  have  examined  solutions  to  the 
strain/relaxation  problem  (sections  4  and  5)  using  both  paths;  they  differ  only  modestly  at  larger  values  of 
Ro,  and  their  qualitative  behaviors  are  the  same.  The  potential-vorticity  path  has  an  additional  complexity 
in  the  strain-impulse  limit  (section  5),  where  its  solution  is  divergent  as  ftmax  -»•  0  (a  limit  that  is  formally 
inappropriate  for  balanced  dynamics  anyway),  although  it  is  well  behaved  as  long  as  Ro~*  >//max  >  1.  In  this 
report,  for  simplicity  we  do  not  show  results  from  the  potential-vorticity  path.  It  is  a  common  experience 
(McWilliams  2003)  that  alternative  balance  approximations  give  modestly,  but  not  importantly,  different 
answers  in  the  regime  of  intermediate  Ro  where  they  are  most  relevant. 
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where  the  quadratic  advection  terms  due  to  the  fluctuations  are  defined  by 

u'M',  +  '  v'M'i  ~  '  (K)2 
EjV"  =  +  u'jdru'e  +  v'idxu'e  -  V/e 

+  u'edru'i+]-v'edxu'i-]-Vev'i 


(9) 


E1V’’  - 


<Wi  +  -/Av'i  +  \ 

+  u'jdrV',  +  —  r'-3/V^  + 
+  «p3rv'-  +  —  if  3xv(  + 


u’iV'i 


u'V 

i  e 


■ 


UN*  = 


-(/)'■  3,.  (rw')  +  m'3, •</>'•  +  -  3  x(v'-0-) 

+  1  (p'^r{ru'e)  +  n-3,-^  +  ^  3 x(v-<#,)  '  ■ 
+  -<#,3r(ra-)  +  <3r0-  +  -  3r(v>') 

Y  Y 


(10) 


(ii) 


The  fluctuation  nonlinear  terms  in  (8)  are  calculated  by  solving  (6)  for  </>'.  evaluating 
(5)  for  n'-  and  v'-,  substituting  into  (9)  (1 1),  and  azimuthally  averaging. 

For  the  vortex  changes  we  make  a  Gradient-wind  Balance  Equation  (BE)  approxima¬ 
tion  (MGM03)  to  exclude  inertia-gravity  oscillations.  In  (8)  this  is  accomplished 
simply  by  dropping  the  radial  acceleration  term,  3 ,(«,•),  in  the  first  equation;  hence, 
rearrangement  yields 


(ut)  =-r[3,(vi)  +  (SJVT) +  ?(«,)] 

*1 

(v,)=|[3,.(0/)  +  (EiVI'>-2^(vf)].  (12) 

Substituting  (12)  into  the  third  equation  of  (8)  and  multiplying  by  —y2  gives 


Y 


-2  f 


-y2dt(c/>)  =  y2(X  3, 

Y 


-  CENV) 

L  C1 


»V 


2^3 t(ve)  -zz£(ue) 
L  Y  9 


r[3,(EiV")] 
y2  - 

+  y  3  ,(rgh(ue)). 


(13) 


where  the  operator  f  is  defined  by 


—2 

f[g]  =  y3r 


Y 


Q 


(14) 
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The  initial  conditions  are  zero  for  all  bracketed  quantities.  The  boundary  conditions 
for  {(pi)  are  9,(0,)  =  0  at  r—0  and  decay  ~r-4  at  the  outer  boundary  for  {(pi)  due  to 
the  external  strain  (section  3)  or  even  more  steeply  ~  e~r/Ld  after  the  strain  abates. 

Equations  (6)  and  (13)  and  their  supporting  relations  comprise  a  well-posed 
initial-  and  boundary-value  problem  for  the  evolution  of  the  fluctuations  and  vortex 
changes  under  the  influence  of  a  specified  external  flow.  They  have  well-defined 
quasigeostrophic  (QG)  and  barotropic  (BT)  dynamical  limits  obtained  by  Ro  -*  0 
and  Lj  — >  oo. 

To  assess  whether  the  straining  and  relaxation  processes  weaken  or  strengthen  the 
vortex,  we  evaluate  norms  based  on  the  area  integrals,  of  the  total  energy  and  potential 
enstrophy  densities,  viz., 

£swe[y,  u,  4>\=\  ( h{u 2  +  v2)  +  Vswe[v,  u,  tp]  -  \hq2,  (15) 

where  h  —  H+tp/g,  q  —  (f+^)/h,  and  £  —  r_1(9r[rv]  —  9^[m]).  Both  of  these  norms 
are  preserved  during  a  conservative  evolution  for  the  vortex,  although  they  are  not 
conserved  when  the  external  strain  flow  is  included.  To  assess  the  changes  in  the 
vortex,  we  subtract  the  initial  values  of  the  norm  associated  with  the  initial  vortex 
(e.g.,  v,  etc.),  and  we  either  include  the  fluctuation  components  (e.g.,  v-,  etc.)  or  not 
in  evaluating  (15)  to  distinguish  between  changes  to  the  fluctuation  part  of  the  vortex 
perturbation  or  to  only  its  azimuthally  averaged  part.  For  example,  excluding  the 
contributions  from  the  external  strain  flow, 

A £swe(t)  =  Sswebi  =  v  +  e2(Vi)(t),  etc.]  -  £we[v,-  =  v,  etc.]  (16) 

is  a  measure  of  the  energy  change  density  that  has  occurred  in  the  azimuthally  averaged 
vortex  structure,  and 

A£'we(0  =  £swe[v  =  v;  =  v  +  ev\{t)  +  e2(v,)(t),  etc.]  -  £swe[v  =  v,-  =  v  +  e2(v, •>(?),  etc.] 

(17) 

measures  the  energy  density  in  the  fluctuations.  We  denote  the  area  integrals  of  (16)  and 
(17)  by  A Eswe  and  A E'sn,e,  respectively,  and  the  analogous  potential  enstrophy  change 
measures  are  AVswe  and  A  Vswe. 


3.  Problem  posing  and  computational  method 

The  initial  vortex  is  chosen  to  have  a  monotonic  relative  vorticity  profile, 


f(r)  =  £ o<r(r/ro)\  (18) 

and  thus  a  finite  circulation  (i.e.,  area  integral  vorticity).  The  constants,  (o  and  r0,  are 
chosen  such  that  the  associated  swirl  profile,  v(r)  from  (4),  has  a  maximum  wind  speed 
V  at  a  radius  of  L.  For  definiteness  in  our  primary  cases,  we  choose  F=4ms-1  and 
L  =  200km.  Along  with  /=10_4s_I,  these  values  imply  Ro—V/fL  —  0.2, 
although  a  related  quantity,  |£|(r  =  0)//has  a  rather  larger  value  of  0.7.  The  associated 
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geopotential  <p  is  calculated  from  (2)  with  an  outer  boundary  condition  of  0  =  0  at 
r  —  Rmax.  We  choose  Rmax  =  2000  km  after  confirming  in  trial  solutions  of  (6)  and 
(13)  that  this  is  sufficiently  distant  for  its  specific  value  not  to  matter.  The  behavior 
shown  below  is  not  qualitatively  sensitive  to  the  particular  profile  choice  in  (18) 
(except  perhaps  in  the  limiting  case  of  a  step-function  profile  where  critical-layer 
interactions  and  quasi-mode  behaviors  may  become  more  prominent). 

We  choose  the  following  as  a  canonical  external  strain  flow  that  acts  only  over  a 
limited  time  interval,  0  <  /  <  tmax : 

u'e  —  —  ctF(/)rcos[2A],  v'  =  +aF(t)r  sin  [2  a],  (ue)  —  (ve)  —  0,  (19) 

where  a  is  the  strain  rate  and 


F(t )  =  sin2 


JTt 
tmax . 


=  0, 


0  <  t  <  t, 
t  —  tmax 


(20) 


is  its  non-dimensional  time  history.  The  cumulative  straining  magnitude  by 
(19)-(20)  after  tmax  is 


d  taF(t)  —  — ^2.  =  (21) 

and  we  choose  m  -C  1  consistent  with  e  1  in  (1).  This  flow  is  zero  at  the  vortex  center, 
is  non-divergent,  and  has  zero  vorticity.  The  associated  external  geopotential  is 

<t>'e  =  y  F(t)r  sin  [2a].  (&>  =  -  l-  (aFrf,  (22) 


which  is  in  (r,  a)  gradient-wind  balance  with  (19).  Again,  the  results  we  report  are  not 
particularly  sensitive  to  the  choice  for  ( ue ,  ve )  as  long  as  its  spatial  scale  is  much  larger 
than  L. 

As  the  centerpiece  of  this  work,  we  focus  on  four  primary  cases:  CYCLONE 
(with  £o  >  0  in  (18)  and  Ro  =  0.2);  ANTICYCLONE  (with  £0  <  0  in  (18)  and 
Ro  =  0.2);  QG  (deleting  all  (9(Ro)  contributions  in  (6)  and  (13)),  which  has  equivalent 
solutions  for  either  sign  of  £0  though  we  do  show  results  for  a  cyclonic  case  with  Co  >  0; 
and  BT  (the  same  as  QG  except  Ld  — >  oo,  hence  y2  — ►  0,  e.g.,  by  taking  the 
limit  H-*-  oo)f,  which  is  directly  comparable  with  the  solutions  in  BG99.  For  the 
non-barotropic  cases  the  layer  depth  is  chosen  to  be  FI—  50  m,  which  yields  an  Lj 
value  of  220  km,  close  to  the  value  of  L  and  consistent  with  parameter  regime  surveyed 
by  Cushman-Roisin  and  Tang  (1990)  and  Polvani  et  al.  (1994).  The  duration  of  the 
strain  period  is  tmax  —  4500  s  (1.25  h);  thus,  it  is  of  brief  duration  (i.e.,  impulsive) 


fBarotropic  dynamics  has  the  same  velocity  evolution  for  all  values  of  Ro.  However,  its  associated 
geopotential  field  is  in  gradient-wind  balance  and  thus  depends  on  Ro.  We  choose  Ro  -y  0  for  the 
BT  case  for  direct  comparability  to  the  QG  solution  with  finite  Lj  and  so  that  (j>jf  is  interpretable  as  the 
streamfunction. 
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MEAN  ENERGY  NORM  MEAN  ENSTROPHY  NORM 


t  V/L  t  V/L 


Figure  2.  History  of  the  normalized  vortex  change  norms  for  (left)  energy,  A Eme(t),  and  (right) 
potential  enstrophy,  A  Vswe(t),  for  the  primary  cases  with  impulsive  strain:  CYCLONE  =  solid  line; 
ANTICYCLONE  =  long-dash  line;  QG  =  dash-dot  line;  and  BT  =  dotted  line.  Normalizations  are 
in  table  1 .  All  curves  start  at  zero  for  t  =  0,  but  the  strain-weakening  period  is  so  brief 
(i.e.,  tmax  V/L  =  0.014)  on  the  scale  of  the  plot  that  they  appear  to  start  at  their  post-strain,  minimum  values. 


compared  to  the  vortex  advective  turn-around  time,  2nL/V  ~  3.6  days,  the  character¬ 
istic  time  for  VRW  evolution  (MGM03).  We  also  choose  a  rather  weak  external  strain 
rate,  a  —  10~7  s-1,  hence  in  =  2  x  10~4  <£  1.  All  of  the  primary  cases  are  run  for  5  eddy 
turn-around  times  (defined  at  the  radius  of  maximum  swirl  velocity),  by  which  point  the 
VRW  relaxation  process  has  asymptotically  approached  its  long-time  limit  in  the  sense 
that  the  energy  and  enstrophy  norms  are  no  longer  changing  significantly.  In  addition, 
other  cases  (defined  in  sections  6  and  7)  are  investigated  to  explore  sensitivities  in  this 
primary  problem  formulation. 

The  numerical  model  used  to  solve  the  tendency  equations  for  the  geopotential 
is  semi-spectral,  with  Fourier  modes  in  the  azimuthal  coordinate  and  second-order, 
centered  finite  differences  in  the  radial  coordinate.  The  nonlinear  terms  are  computed 
by  explicitly  summing  a  convolution  sum  in  Fourier  space.  Inversion  for  the  geopoten¬ 
tial  tendencies  3 ,</>'•  and  3 ,(</>,)  is  accomplished  with  an  LU-decomposition  of  a 
band-diagonal  matrix  derived  from  the  elliptic  operator  in  equations  (6)  and  (16). 
The  tendencies  in  the  nonlinear  terms  on  the  right-hand  side  are  used  from  the  previous 
timestep.  The  numerical  model  is  time-stepped  with  a  fourth-order  Runge-Kutta 
scheme.  A  time  step  of  100  s  with  a  radial  grid  spacing  of  5  km  falls  within  the  CFL 
computational  stability  threshold. 


4.  Strain-induced  net  weakening  and  anticyclonic  dominance 

The  essential  phenomenon  is  illustrated  by  the  histories  of  the  vortex  change  norms  in 
the  four  primary  cases  (figure  2).  Note  that  these  are  area  integrals  of  the  densities  given 
in  (15)  and  (16).  In  all  cases,  the  azimuthally  averaged  energy  and  potential  enstrophy 
decrease  during  the  straining  period  (0  <  t  <  tmax),  and  then  they  partially  recover 
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Figure  3.  (a)  q(r)  normalized  by  flgH  and  (b)  G(r)  normalized  by  fjL  in  CYCLONE  (solid  line)  and 

ANTICYCLONE  (dashed  line). 


during  the  relaxation  period  (tmax  <  t  <  oo).  These  two  periods  are  analyzed  separately 
in  sections  4  and  5.  Note  that  the  changes  in  the  azimuthally  averaged  vortex  are 
much  larger  in  CYCLONE  than  those  in  ANTICYCLONE,  with  QG  in  between 
and  BT  outside  their  range. 

In  advance  of  the  detailed  analysis,  we  can  give  a  simple  heuristic  rationalization 
for  anticyclonic  dominance.  The  external  strain  field  will  induce  a  displacement  \i  of 
vortex  parcels  away  from  their  equilibrium  position.  This  is  a  conservative  process, 
in  particular  of  potential  vorticity.  Within  the  vortex,  the  dominant  potential  vorticity 
is  that  of  the  mean  vortex  itself,  q(r).  Hence  we  expect  potential  vorticity  fluctuations 
that  can  be  estimated  for  small  displacements  as 


Because  of  the  combined  nonlinearities  in  gradient-wind  balance  and  potential  vorticity 
for  finite  Ro,  d,xi  is  very  much  larger  for  cyclones  than  anticyclones  (figure  3a). 
Hence,  for  equal  parcel  displacements  within  the  vortex,  the  q'  fluctuations  will  be 
larger  in  cyclones.  Hence  the  vortex  changes  driven  by  the  fluctuations  will  be  larger. 
Since  the  net  change  is  weakening  (figure  2),  cyclones  will  be  weakened  more,  and 
this  provides  an  evolutionary  path  towards  relatively  stronger  anticyclones.  Note  that 
this  argument  makes  no  specific  reference  to  the  governing  equations  or  the  vortex 
and  external  flow  shapes;  so  we  can  further  expect  that  the  behavior  analyzed  here 
will  occur  under  many  different  circumstances. 


5.  Vortex  weakening  by  impulsive  strain 

To  understand  the  mechanism  of  vortex  weakening  caused  by  an  external  strain  flow, 
we  solve  an  approximate  problem  for  the  vortex  evolution  when  the  straining  event 
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has  a  short  duration  tmax  compared  to  the  free  evolution  of  VRWs  that  occurs  on  the 
advective  time  scale  L/V.  The  size  of  the  net  effect  from  (19)  is  measured  by 

2  tn  f* 

M(t)  =  —  /  F{f) df.  (23) 

‘■max  JO 

M  is  zero  at  /  =  0,  increases  monotonically  for  all  t  >  0,  and  is  equal  to  m  defined  in  (21) 
for  all  t  >  tmax ■  Consistent  with  the  approximation  e  1  (section  2),  we  assume 
that  0  <  m  <£  1,  which  implies  a  small-amplitude  vortex  perturbation  by  the 
external  strain  flow.  The  strain-impulse  approximation  is  a  finite-time  version  of  the 
delta-function  limit  of  impulsive  strain  (i.e.,  tmax  — >  0,  a  -»  oo,  and  m  finite),  and 
its  solution  is  the  same  for  all  times  L/  V '»  t  >  tmax .  More  generally  the  external 
strain  flow  can  have  its  duration  time  tmax  be  either  comparable  to  or  longer  than 
the  vortex  Rossby-wave  time  (cf  sections  6  and  7).  The  processes  of  strain-induced 
weakening  and  wave-induced  recovery  (section  6)  will  then  occur  simultaneously 
and  with  mutual  interaction.  Nevertheless,  the  most  useful  paradigm  for  strain- 
induced  weakening  is  the  response  of  a  stationary  vortex  to  an  impulsive  external 
strain  flow. 

The  impulsive-strain  approximation  allows  us  to  neglect  the  VRW  propagation  terms 
proportional  to  £2  and  G  on  the  left  side  of  (6).  With  a  plane-strain  external  flow  (19), 
the  fluctuation  equation  is 

[v2-y2]^  =  M(/)cos[2k] 

with  S  —  rd,[V/r]  the  strain  flow  for  the  mean  vortex  and  M  the  time  derivative  of  M(t). 
The  solution  of  (24)  has  a  separable  form, 


ij%d_ 
r  dr 


rS 

.  7  . 


-if vr  - 
r  —  hi 

v  \ 


I  1 


(24) 


(j)\r,  X,  t )  —  M(t )  cos[2k]  e//(r),  (25) 

with  (j)\r )  determined  from  (24)  after  substituting  (25)  and  factoring  out  the  ( X ,  t) 
dependences  (i.e.,  the  resulting  right-side  forcing  is  the  term  within  the  outer 
brackets).  In  the  more  easily  interpreted  QG  limit  (Ro  — >-  0),  the  <p\r)  equation 
simplifies  to 


'ld_ 

r  dr 


Hli^f 


(26) 


The  right-side  forcing  in  (26)  vanishes  as  r  ->  0,  oo,  and  it  has  a  single  extremum  in 
between.  For  finite  Ro  values  the  more  general  right-side  forcing  in  (24)  has  a  qualita¬ 
tively  similar  shape.  The  strain-induced  fluctuation  (25)  grows  monotonically  with  an 
azimuthal  pattern  in  quadrature  with  the  strain  flow’s  (22).  The  overall  magnitude  of 
0'  is  0(mfVL). 

The  radial  structure  function  for  the  fluctuations  rp'(r)  is  shown  in  figure  4  for  the 
standard  cases  in  section  3.  It  has  the  shape  of  a  single  dominant  extremum 
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0  2  4  6  8  10 

r/L 

Figure  4.  defined  in  (25)  and  normalized  by  VfL ,  for  the  primary  cases  defined  in  section  3: 

CYCLONE  =  solid  line;  ANTICYCLONE  (with  the  sign  reversed)  =  long-dash  line;  QG  =  dash-dot  line; 
BT  =  dot  line.  The  radius  is  normalized  by  L. 


just  beyond  the  radius  of  maximum  vortex  velocity,  consistent  with  the  shape  of 
5  in  (26).  Compared  to  the  QG  solution,  the  CYCLONE  solution  has  a  much 
stronger  fluctuation  amplitude,  and  the  ANTICYCLONE  solution  has  a  much 
weaker  one.  These  tendencies  are  a  simple  consequence  of  the  larger  magnitudes  of 
rj/f,  f//,  and  Lpy2  in  CYCLONE  amplifying  the  right-side  in  (24)  and  their  smaller 
values  diminishing  the  right-side  in  ANTICYCLONE.  Furthermore,  the  two 
groups  of  forcing  terms  on  the  right  side  of  (24)  are  additive  for  cyclones  with  finite 
Ro,  and  they  have  opposing  signs  for  anticyclones;  this  further  accentuates  the 
reduction  in  the  net  forcing  magnitude.  These  forcing  differences  imply  that  external 
straining  disrupts  cyclones  more  substantially  than  anticyclones  as  Ro  increases, 
through  the  advective  interaction  of  the  external  strain  flow  and  the  mean  vortex. 
It  is  interesting  that  0'(/-)  (hence  also  qah)  changes  sign  in  the  outer  region  for 
ANTICYCLONE;  this  will  prove  to  be  significant  in  the  response  to  sustained 
strain  (section  8).  Finally,  <p\r)  is  modestly  larger  in  the  BT  solution,  compared 
to  QG;  both  solutions  have  the  same  right-side  forcing  in  (26),  but  BT  has  a  larger 
fluctuation  response  (i.e.,  cp/q  ratio)  than  QG  since  it  lacks  the  vortex-stretching  term 
oc  Lp2  in  q. 

The  vortex  change  equation  (13)  has  a  composite  separable  solution  form,  viz., 

mr,t)  -  M\t )  (&M(r)  +  (0)M(r).  (27) 

Its  first  term  is  0(nrfVL),  and  it  exhibits  monotonic  growth  until  t  —  tmax  after  which  it 
is  steady.  Its  second  term  is  0{m2fVL/(ftmax)2).  It  is  at  least  as  large  as  the  first  one 
in  the  impulsive-strain  limit,  but  its  temporal  structure  is  a  single-signed  pulse  in  the 
interval  between  vanishing  at  t  —  0  and  t  >  tmax.  Thus,  the  second  term  represents  a 
fully  reversible  response  during  the  straining  event,  and  the  persistent  vortex  change 
caused  by  straining  is  entirely  represented  by  the  first  term,  ( <p)M . 
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Consistent  with  the  solution  form  (27),  we  decompose  the  fluctuation  nonlinear  terms 
in  (13)  as  follows: 
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The  resulting  radial  boundary  value  problems  for  the  components  of  (0)  in  (27)  are 


hL(_7i-i 
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(29) 


The  velocities  for  the  vortex  change  (12)  have  an  analogous  decomposition  into 
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The  resulting  relations  for  their  radial  functions  are 


(u)j  = 


~  |  (/(%  +  ), 


i  (mi 
l\  d'- 


(31) 


for  j  —  M,M  and  Sjj  the  discrete  delta  function  (i.e.,  5=1  for  equal  indices  and  5  =  0 
otherwise).  The  vortex  change  velocity  magnitudes  are  0(m2  V/(  ftmax))  and 
0(m2  V/(  ftmax)3),  respectively,  for  the  two  (u)j  components  in  (30),  and  they  are 
0(m2V)  and  0(m2V/(ftmax)2),  respectively,  for  the  two  (v)j  components.  They  are  all 
small  compared  to  v  by  the  factor  nr,  and  they  are  enhanced  in  magnitude 
during  the  strain-impulse  phase  by  various  factors  of  1  /(ftmax).  After  the  strain  impulse 
has  ended  (i.e.,  for  t  >  tmax),  the  only  nonzero  velocity  component  is  (v)M  with  a 
magnitude  of  0(m2V). 

In  the  QG  case,  A =  Af’Jy  —  0,  and  NUM,  although  nonzero,  does  not  contribute  to 
the  leading-order  relations  for  ( u ),  and  (v).  The  limiting  form  for  (29)  is 


"Id 
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with  the  eddy  momentum  flux  forcing, 
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(33) 


The  associated  QG  velocity  components  are 


(u)  =  0, 


(v> 


1  d(0)M 

/  dr 


(34) 


The  radial  structures  that  survive  at  the  end  of  the  strain-impulse  period,  (<p)M  and 
{v)M,  are  shown  in  figure  5  for  the  standard  cases  in  section  3.  Their  profiles  are 
generally  similar  among  the  cases.  (</>)// is  positive  within  the  vortex  core,  opposite  to 
the  sign  of  <p/f.  (v)  is  also  opposite  in  sign  to  v  near  its  extremum  at  r  =  L.  These 
indicate  a  weakening  of  the  vortex  during  the  straining  event.  The  magnitudes  of  the 
vortex  change  profiles  are  larger  for  CYCLONE  than  ANTICYCLONE  and  smaller 
for  QG  than  for  BT. 

The  dominant  forcing  term  for  (</>)M  in  (29)  is  A fy,  and  its  largest  components 
are  those  analogous  to  the  QG  ones  in  (33).  The  dominance  of  these  terms  also 
occurs  for  (v)  in  (31).  The  principal  source  of  cyclone-anticyclone  differences  is  the 
G(Ro)  forcing  terms  involving  the  advective  interaction  of  the  external  strain 
flow,  the  mean  vortex,  and  the  strain-induced  fluctuation,  <//  (i.e.,  not  those  involving 
quadratic  products  of  the  fluctuations  nor  those  of  0{Ro2)).  For  example,  a  member 
of  the  important  O(Ro)  class  of  terms  in  A fr  in  (28)  is  —S^'/rj^r.  In  the  case  of  cyclones 
these  extra  terms  generally  reinforce  the  dominant  QG  forcing  terms,  whereas  for 
anticyclones  they  have  opposing  signs  and  reduce  the  net  right-side  forcing. 
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Figure  5.  (a)  (<A)mO')  normalized  by  VfL  and  (b)  [v)M(r)  normalized  by  V  for  the  standard  cases  defined 

in  section  3:  CYCLONE  =  solid  line;  ANTICYCLONE  (with  the  sign  reversed)  =  long-dash  line;  QG  = 
dash-dot  line;  BT  =  dot  line.  The  radius  is  normalized  by  L. 


The  BT  changes  in  figure  5  are  larger  than  the  QG  changes  (with  finite  Ld )  for  the 
same  reason  as  the  difference  in  0'  amplitudes,  viz.,  the  elliptic  operator 
inversion  in  (29)  gives  a  larger  response  for  similar  right-side  forcing  without  the 
second  term  associated  with  vortex  stretching.  The  resulting  difference  is  larger 
for  the  vortex  changes  than  the  fluctuations  since  the  <p/q  ratio  is  larger  for  the 
vortex-change  response  because  the  left-side  operator  lacks  the  terms  associated  with 
the  azimuthal  second  derivatives  in  (24). 

We  use  several  norms  to  assess  the  vortex  evolution  under  strain  (section  2).  The  fluc¬ 
tuation  and  bracket  norms  from  (15)  are  positive  quadratic  functionals  of  <//  and  (0), 
respectively.  Consequently,  these  norms  exhibit  a  combination  of  both  transient  spiking 
and  monotonic  growth  during  the  strain-impulse  period,  reflecting  the  pattern  of 
perturbation  growth.  To  assess  the  overall  vortex  change,  however,  we  use  the 
full  SWE  energy  and  potential  enstrophy  norms.  The  SWE  energy  change  for 
the  azimuthal-averaged  vortex  (16)  is  evaluated  for  these  solutions.  We  focus  on  the 
irreversible  change  that  occurs  for  t  >  tmax,  and  we  make  the  approximation  m  <<:.  1. 
The  result  is 

A Eswe  ^2tv  j“  drr(hv{v)  +  ^ v2)(0) (35) 

evaluated  at  t  —  tmax.  The  magnitude  of  A Eswe  is  0{m2HL2V 2).  The  third  term  in  the 
integrand  is  smaller  than  the  first  two  terms  by  O(Ro).  Thus,  the  essential  cause  for 
vortex  weakening  by  external  strain,  as  measured  by  this  energy  norm,  is  contained 
within  the  first  radial  boundary-value  problem  in  (29)  and  its  associated  velocity 
relations  (31)  that  yield  radially  integrated  negative  correlations  between  the  mean 
vortex  and  its  strain-induced  change  for  both  the  mass  field  (</>(</>}  <  0  on  average) 
and  the  azimuthal  velocity  (v(v)  <  0).  The  values  of  A Eswe  for  the  standard  cases  in 
section  3  are  listed  in  table  1.  In  all  the  cases,  strain  weakens  the  mean  vortex 
energy,  and  it  does  so  much  more  for  a  cyclone  than  an  anticyclone. 

The  analogous  measure  for  the  energy  change  associated  with  the  azimuthal 
fluctuations  is  defined  in  (17).  Again  making  the  approximation  m  1  and  noting 
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Table  1.  Strain-impulse  norms. 


CASE 

A  Eswe 

^ Kswe 

AKwe 

CYCLONE 

-3.50 

1.83 

-54.2 

16.7 

ANTICYCLONE 

-1.04 

0.35 

-1.74 

0.48 

QG 

-3.02 

1.21 

-8.84 

2.32 

BT 

-4.41 

1.50 

-5.03 

2.10 

The  energy  and  potential  enstrophy  norms  evaluated  from  (35) — (47)  at  /  =  tmax  and  normalized  by  lirHLrV 2  for  energy  and 
mV2 jg  H  for  enstrophy. 


that  all  odd  powers  of  the  azimuthal  fluctuations  in  Eswe  have  zero  azimuthal  average, 
the  appropriate  norm  is 

AKwe  ~  j  dx  Jo  drr(^hi(u,2  +  v'2)  +  ^0,2'j  (36) 

evaluated  at  t  —  tmax.  The  magnitude  of  A Eswe  is  0(m2HL2V2),  i.e.,  the  same  as  A Eswe. 
The  final  term  in  the  integrand  is  smaller  than  the  first  three  terms  by  O(Ro).  Since  the 
dominant  first  three  terms  are  positive-definite  in  A E!swe,  a  strain-impulse  acts  primarily 
to  increase  the  energy  associated  with  azimuthal  fluctuations.  The  values  of  A E'swe  for 
the  standard  cases  in  section  3  (table  1)  are  all  positive,  indicating  the  generation  of 
fluctuations  by  straining.  The  fluctuation  energy  change  is  much  greater  for  cyclones 
than  for  anticyclones. 

The  SWE  potential  enstrophy  change  norms  (15)  at  t  —  tmax,  with  the  same 
approximations  as  for  the  energy  norms,  are  the  following: 


A  Vswe  ^  2JT 


aV... 
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p2  tc  noc 

/  / 

Jo  Jo 


dr%2{,) 

2 g 


g2H 


<?(</>» 


/  gH)\f  gH 


dr—U'2+f<P'2-2c,W). 
2  g2h 


(37) 


The  magnitude  for  both  these  norms  is  0(m2V2/g2H).  The  arrow  in  (37)  indicates  the 
leading-order  expression  in  the  QG  limit  for  the  enstrophy  vortex  change  norm;  its 
magnitude  is  formally  smaller  than  the  magnitude  of  the  first  expression  by  a  factor 
of  Ro,  since  evaluating  the  first  expression  with  q  —f/H  gives  a  trivial  result  with  {<j>) 
and  (f)  obtained  from  (32)-(34). 

The  vortex  change  and  asymmetric  fluctuation  enstrophy  norms  in  table  1  indicate 
the  same  relative  differences  between  cyclones  and  anticyclones.  The  enstrophy 
differences  are  even  larger  than  for  energy,  mainly  because  the  h  weights  in  (37) 
greatly  accentuate  them.  The  BT  energy  changes  are  larger  than  the  QG  ones, 
and  the  enstrophy  changes  are  somewhat  smaller.  These  outcomes  are  based  on  the 
competing  effects  of  larger  barotropic  amplitude  (figures  4  and  5)  —  favoring  larger 
BT  changes  —  and  the  absence  of  the  effects  of  changes  in  depth  (i.e.,  terms 
involving  <p)  in  the  barotropic  form  of  (35) — (37),  whose  signs  in  QG  act  to  reinforce 
the  remaining  terms  —  favoring  larger  QG  changes. 

In  summary,  an  external  strain  flow  acts  to  weaken  a  vortex,  and  this  is  most 
cleanly  demonstrated  with  the  strain-impulse  approximation  used  in  this  section. 
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External  strain  induces  larger  vortex  changes  for  cyclones  than  anticyclones  and  larger 
changes  for  barotropic  vortices  than  quasigeostrophic  ones.  The  mean  energy  and 
enstrophy  change  norms  in  table  1  summarize  the  net  vortex  weakening. 


6.  Vortex  recovery  by  VRWs 

This  section  describes  how  a  vortex  recovers  from  a  strain-induced  perturbation, 
and  it  quantifies  the  net  effect  of  a  straining  event  on  the  vortex.  The  straining 
period  specified  in  section  3  has  a  small  duration  tmax  compared  to  the  characteristic 
advective  time  scale  L/V  for  the  vortex  relaxation  process.  The  imposed  strain  flow 
induces  an  azimuthal  wavenumber-two  fluctuation  at  t  —  tmax  whose  pattern  is  separ¬ 
able  in  r  and  X  as  in  (25).  It  has  the  same  sign  for  all  r  values  except  for  the  radial 
sign  change  for  ANTICYCLONE  that  occurs  outside  the  vortex  core  (at  r  «  2.5  in 
figure  4).  Thus,  we  can  say  that  the  initial  value  of  a  ‘local’  radial  wavenumber  on 
scales  smaller  than  the  vortex  size  is  zero,  i.e.,  k(0)  —  0.  This  pattern  can  be  seen  in 
the  upper  left  panel  of  figure  6  for  qah  defined  in  (7).  During  the  relaxation  phase, 
tmax  <  t  <  oo,  the  mean  angular  velocity  profile  Q(r)  causes  the  familiar  differential 
rotation  around  the  vortex  whereby  the  spiral  arms  wrap  up,  propagate  outward, 
and  ultimately  filament  (Montgomery  and  Kallenbach  1997,  MGM03);  this  is  shown 
for  CYCLONE  in  the  subsequent  panels  of  figure  6.  This  behavior  is  represented  by 
a  k(t)  that  grows  monotonically  in  time  in  association  with  an  increasing  number  of 
sign  changes  in  qah  along  all  radial  lines,  and  it  is  also  associated  with  systematically 
decaying  fluctuation  geopotential  and  velocity  amplitudes. 

Figure  7  shows  \qah\{r,t)  for  CYCLONE  and  ANTICYCLONE.  The  patterns 
indicate  an  outward  radial  propagation  of  the  asymmetric  fluctuations  in  both  cases. 
The  stagnation  radius  (i.e.,  the  limit  of  outward  propagation)  for  the  primary  wave 
in  CYCLONE  is  at  approximately  1.8L.  This  is  noticeably  larger  than  the  stagnation 
radius  for  the  ANTICYCLONE,  which  is  1.3L,  and  also  larger  than  in  QG  and  BT, 
which  both  have  a  stagnation  radius  of  1.571.  The  stagnation  radius  is  controlled  by 
the  radial  group  velocity  of  the  VRWs.  Within  the  context  of  waves  in  slowly  varying 
media  (i.e.,  a  WKB  approximation),  the  radial  group  velocity  (equation  (62)  of 
MGM03)  is  controlled  by  the  mean-vortex  potential  vorticity  gradient  G(r)  and 
the  radial  wavenumber  k.  In  CYCLONE,  G  is  negative  everywhere  since  q(r)  is  mono¬ 
tonically  decreasing  (figure  3b),  and  k  becomes  positive.  Since  the  radial  group  velocity 
is  negatively  correlated  with  the  product  of  G  and  k.  this  yields  a  positive  value,  and 
thus  an  outwardly  directed  propagation.  In  ANTICYCLONE,  G  and  k  both  have 
opposite  signs  compared  to  CYCLONE,  but  their  product  is  the  same,  so  this  case 
also  exhibits  an  outward  group  velocity. 

From  (4)  and  (7),  G  can  be  expanded  into 


(38) 


When  terms  of  C9(Ro)  are  quantitatively  important,  (f,  fj)  are  larger  for  cyclones 
because  their  component  terms  are  additive,  whereas  they  are  smaller  for  anticyclones 
because  the  terms  have  opposite  signs  and  thus  partially  cancel.  In  contrast  h  is  smaller 
for  cyclones  for  an  analogous  reason,  and  this  change  is  dominant  compared  to  fj  in 
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Figure  6.  Swirl  plots  for  qah(x,y)  normalized  by  m/K/L  in  CYCLONE  at  several  times  tV/L  during  the 
relaxation  phase  after  impulsive  strain. 
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Figure  7.  Hovmiiller  plots  of  |]  (r,  /)  normalized  by  mfV/L  after  impulsive  strain  in  (a)  CYCLONE 

and  (b)  ANTICYCLONE.  The  contour  interval  is  0.5  in  (a)  and  0.1  in  (b). 


making  q  larger:  9,</>  is  larger  for  cyclones  because  of  gradient-wind  balance  (2),  and  for 
both  vortex  parities  it  combines  additively  with  9,.£  with  a  further  weighting  by  q  in  (38). 
As  a  result  of  these  differences,  the  magnitude  of  G  is  significantly  larger  for  cyclones 
than  anticyclones,  so  are  the  group  velocity  and  the  stagnation  radius  (figure  3b). 
The  outward  radial  group  velocity  causes  a  rapid  evacuation  of  fluctuation  potential 
vorticity  in  the  core  of  the  vortex  that  occurs  in  conjunction  with  the  spiral  wind-up 
of  the  primary  wave  packet.  The  core  evacuation  is  also  due  to  the  outwardly  directed 
group  velocity,  and  it  begins  fairly  early.  The  outward  propagation  of  asymmetric 
potential  vorticity  is  evident  in  the  spiral  plots  of  qah{r,  X)  in  the  core  of  both  cyclones 
(figure  6)  and  anticyclones. 

Because  of  the  conservation  of  wave  action  and  wave  activity  for  VRWs  (MGM03), 
amplitude  decay  of  the  primary  wave  packet  is  expected  for  an  outwardly  propagating 
VRW.  Comparing  the  net  change  in  qah  during  the  relaxation  phase,  there  is  clearly 
a  greater  amount  of  fluctuation  decay  in  CYCLONE  than  in  ANTICYCLONE. 
The  magnitude  of  the  strain-induced  fluctuation  \<fi'i\(tmax)  in  CYCLONE  is  larger 
than  in  ANTICYCLONE,  with  the  QG  and  BT  cases  lying  in  between  (figure  4). 
The  \qab\(tmax)  differences  among  the  cases  are  even  larger,  with  a  ratio  of  about 
a  factor  of  5  between  CYCLONE  and  ANTICYCLONE.  In  the  absence  of  any 
dissipation  for  the  fluctuations,  the  combination  of  greater  cyclonic  qab  after  straining 
and  its  greater  relative  decay  during  the  relaxation  phase  suggests  there  must  be  a  much 
greater  strengthening  of  the  cyclonic  vortex  through  eddy-mean  interaction  in  (13). 
Since  the  fluctuation  is  originally  created  by  transfer  from  the  mean  vortex 
(i.e.,  strain  weakening),  this  VRW  relaxation  strengthening  is  an  opposite  effect.  As 
evident  in  figure  2,  however,  the  strengthening  is  less  than  the  weakening  in  all  cases. 

Figure  8  shows  the  temporal  decay  of  the  fluctuation  geopotential  amplitude  for  the 
four  principal  cases  at  the  radius  of  maximum  mean  potential-vorticity  gradient  after 
the  external  strain  is  turned  off.  The  evolution  of  the  geopotential  at  the  radius  of 
maximum  mean  potential-vorticity  gradient  serves  to  identify  whether  the  relaxation 
response  is  dominated  by  a  single  VRW  quasi-mode  (with  exponential  decay; 
Schecter  et  al.  2002)  or  a  collection  of  sheared  (continuum)  VRWs  that  decay  algebrai¬ 
cally  with  time  (MGM03).  From  figure  8  it  is  evident  that  the  decay  is  only  persistently 
exponential  in  BT.  Decay  in  QG  is  less  persistently  exponential,  if  at  all. 
CYCLONE  and  ANTICYCLONE  exhibit  highly  non-exponential  decay.  These  results 
suggest  that  the  sheared  VRW  paradigm  is  the  most  appropriate  one  for  the  finite-Ro 
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Figure  8.  Decay  of  \<j>\  at  the  radius  of  maximum  mean  potential  vorticity  gradient  (r  rs  0.6L)  during  the 
vortex  relaxation  phase  (starting  at  t=  tmax),  normalized  by  its  value  at  tmax.  CYCLONE  =  solid  line; 
ANTICYCLONE  =  long-dash  line;  QG  =  dash-dot  line;  BT  =  dotted  line. 


Figure  9.  Vortex  changes  during  the  VRW  relaxation  phase  following  impulsive  straining  (i.e.,  between 
t  =  31.4  and  t  =  tmax ):  (a)  (0)(r)  normalized  by  m2VfL  and  (b)(v)(r)  normalized  by  nrV.  CYCLONE  =  solid 
line;  ANTICYCLONE  (with  the  sign  reversed)  =  long-dash  line;  QG  =  dash-dot  line;  BT  =  dotted  line. 


cyclonic/anticyclonic  vortex  recovery  phase.  In  general,  we  expect  that  the  relaxation 
response  will  involve  a  superposition  of  VRWs  and/or  quasi-modes,  but  we  have  not 
made  any  attempt  to  decompose  the  fluctuation  field  into  individual  ‘pure  tones’. 

Figure  9  shows  the  net  geopotential  vortex  change  (</>}  at  /  =  31.4L/Fminus  its  value 
at  the  end  of  the  straining  period,  t  =  tmax.  A  bimodal  structure  in  {<j>)  is  clearly  evident 
for  all  four  cases,  with  a  large  negative  value  inside  the  radius  of  maximum  wind  for  v(r) 
and  a  smaller  positive  value  outside.  The  negative  {(/>)  corresponds  to  a  greater  mean 
central  depression  in  h  for  a  cyclone  and,  by  gradient-wind  balance,  a  stronger 
vortex.  Since  the  cyclonic  {</>}  exhibits  a  peak  amplitude  that  is  approximately  twice 
that  of  the  anticyclonic  case  (plotted  here  with  a  sign  reversal),  this  implies  that  the 
cyclone  undergoes  a  greater  recovery  from  the  strain  weakening.  Note  that  the  peaks 
in  {(/>)  for  the  cyclone  are  slightly  outside  the  others  due  to  the  larger  radial  propagation 
of  the  corresponding  VRWs.  Again,  the  QG  and  BT  results  typically  lie  in  between 
CYCLONE  and  ANTICYCLONE. 

From  equation  (33)  of  MGM03,  relating  the  time  derivative  of  balanced  fluctuations 
to  the  wave-induced  acceleration  of  the  azimuthal  velocity,  the  larger  response  in  (v)  in 
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the  cyclone  seen  in  figure  9  is  traceable  to  the  larger  radial  propagation  in  the  cyclone  in 
an  approximate  pseudo-momentum  rule  in  the  limit  of  small  but  finite  Ro.  For  the  BT 
case  with  non-divergent  flow,  the  physical  interpretation  of  this  wave-induced  accelera¬ 
tion  is  a  familiar  one.  As  the  mean  angular  velocity  causes  the  fluctuation  to  wrap 
around  the  vortex,  the  phase  lines  tilt  in  the  direction  of  the  shear.  This  generates 
an  eddy  angular  momentum  flux  that  is  convergent  on  the  inward  side  of  the  waves, 
and  divergent  on  the  outward  side.  Since  the  mean  tangential  velocity  is  negatively 
correlated  with  the  divergence  of  the  eddy  momentum  flux,  this  implies  an  increase 
in  (v)  on  the  inward  side  of  the  waves  and  a  decrease  on  the  outward  side.  This  feature 
is  seen  in  figure  9  for  all  the  four  cases,  and  implies  a  vortex  strengthening  in  the  core. 
Note  that  the  positive  inner  peak  is  larger  than  the  negative  one  outside  the  core. 
Since  the  total  angular  momentum  of  the  fluid  must  be  conserved  during  the  relaxation 
phase,  the  VRWs  only  redistribute  angular  momentum,  r(v).  Thus  at  the  smaller  inner 
radius,  the  vortex  change  must  be  larger  in  magnitude  to  balance  the  opposite-sign 
change  at  the  larger  outer  radius.  Since  the  majority  of  the  wave-mean  flow  interaction 
occurs  near  the  stagnation  radius,  the  zero  value  in  (v)  will  occur  near  this  location. 
The  profile  of  (v)  for  CYCLONE,  with  its  larger  VRW  stagnation  radius,  is  thus  shifted 
outwards  compared  to  the  other  cases. 

The  energy  and  potential  enstrophy  vortex  change  norms  ( 1 5)  ( 1 6)  quantify  the  area- 
integrated  vortex  change  (figure  2).  These  measures  are  approximately  controlled  by  (<p) 
and  (v)  as  indicated  in  (35)  and  (37).  After  the  initial  (and  here  brief)  strain-weakening 
in  these  norms,  both  energy  and  enstrophy  plots  increase  rather  sharply  at  early  times  as 
the  recovery  begins,  but  asymptote  after  a  time  of  about  t  —  8 L/V,  corresponding  to  a 
little  more  than  one  eddy  turn-around  time.  After  this  point,  the  fluctuations  have 
become  significantly  wrapped  around  the  vortex  so  that  their  phase  lines  are  almost 
aligned  with  the  mean  swirl  flow,  so  the  eddy  momentum  fluxes  become  small.  The 
recovery  phase  for  CYCLONE  can  be  seen  to  be  larger  in  magnitude  for  both  energy 
and  enstrophy  compared  to  ANTICYCLONE,  as  can  be  anticipated  by  the  larger 
peaks  in  {(p)  and  (v)  during  this  time.  Since  the  strain-induced  weakening  is  so  much 
larger  for  CYCLONE  despite  this  greater  strengthening,  its  late- time  asymptotic 
values  are  still  more  negative  than  for  ANTICYCLONE  (with  the  QG  case  in  between), 
which  indicates  a  greater  net  weakening  for  CYCLONE. 

The  net  vortex  changes  in  (</>)  and  (v)  over  both  the  straining  period  and  relaxation 
phase  are  in  figure  10.  Note  that  (tp)  shows  larger  positive  values  for  CYCLONE  than 
ANTICYCLONE,  and  in  fact  for  CYCLONE  it  is  almost  exclusively  positive  which 
demonstrates  weakening  of  the  vortex  at  almost  all  radii.  In  (v}(r),  ANTICYCLONE 
has  larger-magnitude  peaks  than  CYCLONE,  which  would  suggest  more  anticyclonic 
weakening.  However,  the  mean  energy  and  enstrophy  changes  depend  more  heavily 
on  the  changes  in  (</>)(r)  where  the  cyclone  is  more  heavily  weakened.  The  large  peaks 
in  both  (<p)  and  (v)  for  the  BT  case  reflect  the  difference  in  the  (p/q  strain-response 
ratio  due  to  the  absence  of  vortex-stretching  term  when  y  —  0  (section  5).  The 
manifestation  of  the  larger  vortex  change  profiles  in  the  energy  and  enstrophy  change 
norms  (i.e.,  larger  A Eswe,  smaller  AVswe;  figure  2)  reflects  a  subtle  competition  between 
the  velocity  and  geopotential  components  in  (35)  and  (37). 

In  summary,  during  the  relaxation  phase,  the  perturbed  vortex  exhibits  the  following 
behaviors.  For  a  fluctuation  in  the  absence  of  further  external  straining,  there  is  radial 
propagation;  azimuthal  wrap-up  and  (r,  X)  phase-line  tilting  down  the  gradient  of  Q(r); 
and  rapid  amplitude  decay  in  fluctuation  geopotential  and  velocity.  The  effect  of  the 
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r/L  r/L 

Figure  10.  Net  vortex  change  in  (a)  (</>)(/■)  normalized  by  m2  VfL  and  (b)  (v)(r)  normalized  by  m2  V  over  the 
full  time  integration,  including  both  the  impulsive  straining  and  relaxation  phases.  CYCLONE  =  solid  line; 
ANTICYCLONE  (with  the  sign  reversed)  =  long-dash  line;  QG  =  dash-dot  line;  BT  =  dot  line. 


mean  eddy  fluxes  on  the  vortex  is  to  strengthen  it  in  both  energy  and  enstrophy 
during  this  period  of  strong  phase-line  rotation.  Based  on  an  extensive  survey  of 
vortex-relaxation  solutions  (not  reported  here),  we  have  found  no  pervasive  differences 
between  the  relative  efficiency  of  cyclonic  and  anticyclonic  strengthening  due  to 
this  process.  This  is  consistent  with  the  similar  proportional  changes  seen  here  in 
CYCLONE  and  ANTICYCLONE.  Thus,  we  conclude  that  the  essential  cause  of 
anticyclonic  dominance  is  the  smaller  strain  weakening  for  anticyclones  (section  5). 
All  of  these  fluctuation  and  vortex-change  behaviors  are  consistent  with  previous 
studies  of  VRWs  and  quasi-modes  (e.g.,  MGM03;  Schecter  et  al.  2000). 


7.  Sensitivities 

7.1.  Deformation  radius 

For  the  primary  cases  presented  in  section  3,  Ro  =  0.2,  and  Ld/L  —  1.11.  Elere  we 
examine  cases  with  a  Ro  that  is  half  as  large,  Ro  =  0.1,  to  allow  for  smaller  Ld 
values  that  would  otherwise  bottom-out;  i.e.,  the  cyclonic  h  would  vanish  in  the 
vortex  core.  Specifically,  we  compare  cyclonic  and  anticyclonic  solutions  for  three 
Ld/L  values  equal  to  0.78,  1.11,  and  1.57  (corresponding  to  mean  depths  of  H  —  25, 
50,  and  100  m).  When  LdjL  is  large  enough,  the  flow  becomes  essentially  barotropic, 
a  regime  which  has  dynamical  symmetry  between  cyclones  and 
anticyclones.  We  would  expect  then  that  as  we  move  further  from  the  barotropic 
regime  (i.e.,  decreasing  Ld/L),  we  should  begin  to  see  a  clear  distinction  between 
cyclones  and  anticyclones. 

Figure  1 1  shows  (v)(r),  for  the  three  LdjL  values.  As  LdjL  decreases,  not  only  do  both 
cyclones  and  anticyclones  have  a  reduced  peak  value,  but  the  cyclonic  and 
anticyclonic  profiles  diverge  as  expected.  This  result  is  reflected  in  the  mean  enstrophy 
norm  (i.e.,  ( 1 5)  ( 1 6)  and  (37))  and  summarized  in  table  2  for  the  ratio  of  the  late-time 
(/  =  31.4L/L)  enstrophy  change  norms,  A  Vswe,  for  cyclonic  and  anticyclonic  mean 
vortices.  Since  all  cases  demonstrate  net  vortex  weakening,  their  mean  norms  are 
negative  (with  larger  negative  values  for  larger  weakening).  Thus  ratios  greater  than 
one  indicate  anticyclonic  dominance,  which  increases  as  the  ratio  grows.  We  have 
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Ld/L  =  1-57 


Ld/L  =  1.11 


Ld/L  =  0.78 


Figure  11.  Changes  in  (v)(r)  normalized  by  m2V  at  i  =  31.4L/K  for  CYCLONE  (solid  line)  and 
ANTICYCLONE  (long-dash  line,  with  the  sign  reversed)  with  impulsive  strain  forcing  and  different  values 
of  Ld/L. 


Table  2.  Cyclone/anticyclone  norm  ratios. 


Ro 

Ld/L 

A  Vswe  ratio 

0.1 

1.57 

1.86 

0.1 

1.11 

3.29 

0.1 

0.78 

10.3 

0.2 

1.11 

20.9 

The  ratio  of  late-time  (i.e.,  t=3lAV/L)  cyclonic  and  anticyclonic  vortex  change 
enstrophy  norms  with  impulsive  strain  and  different  values  of  L^/L. 


also  included  the  value  for  the  primary  cases  given  in  sections  4  and  5  for  comparison 
(where  Ro  =  0.2). 

In  the  norms  the  term  involving  (v>  has  the  most  sensitivity  to  changes  in  Llh 
although  the  vortex  change  geopotential,  p ),  also  has  an  important  contribution. 
The  energy  norm  is  not  so  monotonic  in  its  sensitivity  to  changes  in  Ld/L, 
so  it  is  thus  a  less  useful  integral  measure  to  compare  cyclonic  and  anticyclonic 
dependences  on  L^/E. 

7.2.  Rossby  number 

Our  investigation  is  limited  in  Ro  by  the  disappearance  of  cyclone-anticyclone 
differences  as  Ro  -*  0  (i.e.,  QG)  and  the  vanishing  of  h(r  —  0)  as  Ro  increases  in 
cyclones  for  a  given  value  of  L^L.  In  the  SWE  with  the  mean  vortex  profile  (18)  and 
Ld/L  ~  1,  this  upper  bound  on  Ro  comes  into  force  before  the  centrifugal 
instability  boundary  (i.e.,  £(r  =  0)  =  — f)  for  anticyclones.  Within  this  Ro  range 
we  have  found  no  surprises.  Typically  CYCLONE  and  ANTICYCLONE  results 
diverge  from  QG  in  opposite  ways  as  Ro  increases  (e.g.,  figure  4  el  seq.),  and  their 
differences  vanish  smoothly  as  Ro  decreases  (e.g.,  table  2). 

7.3.  Strain  duration 

Sections  3-5  focus  on  cases  with  only  a  brief  duration  of  the  external  strain  flow, 
tmax  <£  L/V.  As  also  done  in  BG99,  we  have  examined  various  solutions  with 
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tmax~L/V  and  found  no  qualitative  changes  in  the  phenomena  but  an  increased 
complexity,  because  the  strain-induced  fluctuations  are  continuously  generated  while 
the  VRW  wrap-up  is  in  progress;  furthermore,  the  vortex  change  norms  become 
more  oscillatory  in  time.  However,  there  are  interesting  differences  when  the  external 
strain  is  sustained  in  time,  tmax  L/V,  and  this  regime  is  examined  in  section  8. 

7.4.  Wave  nonlinearity 

The  model  we  are  using  is  a  quasi-linear  one:  the  fluctuation  dynamics  are  linear, 
but  the  mean  vortex  changes  are  due  to  averages  of  nonlinear  fluctuation  products. 
We  wish  to  address  the  extent  to  which  this  linear  wave  approximation  is  appropriate 
for  the  primary  cases  presented  in  this  article,  and  also  to  establish  the  amplitude  of 
the  external  strain  flow  at  which  the  nonlinear  effects  become  important.  To  this  end 
we  add  in  the  nonlinear  wave-wave  interactions,  in  particular  the  interactions  between 
wavenumber-two  fluctuations  to  produce  wavenumber-four  fluctuations,  which  in  turn 
interacts  with  itself  to  contribute  to  the  vortex  change  and  with  the  wavenumber-two 
fluctuation  back  onto  the  wavenumber-two  fluctuation. 

We  examine  the  QG  case  as  the  simplest  and  the  most  straightforward  one.  In  the 
primary  QG  case  (section  3)  but  with  added  wave-wave  terms,  there  are  essentially 
no  important  differences  due  to  wave-wave  nonlinearity.  The  amplitude  of  the  induced 
wavenumber-four  fluctuation  is  about  four  orders  of  magnitude  smaller  than  the 
primary  wavenumber-two  fluctuation,  so  this  result  is  not  surprising.  Consistent  with 
this  result  as  a  guideline,  it  is  not  until  the  external  strain  flow  has  an  amplitude 
a  that  is  about  2000  times  larger  than  that  used  in  the  primary  case  (section  3) 
that  the  wavenumber-four  and  -two  amplitudes  become  comparable,  and  thus 
the  wave-wave  effects  become  more  significant.  Note  that  with  this  a  —  2  x  10-4  s  , 
the  external  strain  velocity  (19)  at  r  —  L  is  10  times  larger  than  the  mean  vortex  V. 
(Alternatively,  with  a  larger  strain  duration  tmax  a  smaller  a  could  achieve  the  same 
effect.)  This  shows  that  the  wave  nonlinearity  is  significant  only  in  extreme  cases. 

Since  the  wavenumber-four  pattern  has  twice  as  many  oscillations  azimuthally 
around  the  vortex,  the  combined  fluctuation  pattern  strengthens  in  some  places 
and  weakens  in  others  compared  to  the  wavenumber-two  pattern  alone.  Also,  the 
wavenumber-four  peak  is  located  further  out  radially  than  the  wavenumber-two 
peak,  so  the  combined  pattern  is  a  modulated  wavenumber-two  pattern.  It  does, 
however,  windup  spirally  as  before,  with  the  underlying  behavior  being  similar. 
Further,  the  mean  norms  display  only  small  differences  from  the  primary  QG  case, 
with  approximately  a  6%  change  in  end-state  behavior  for  the  mean  energy,  and 
about  2%  change  for  mean  enstrophy,  where  both  changes  reflect  less  net  weakening 
of  the  vortex.  Thus  even  at  much  larger  external  strain  flow  amplitudes,  the  wave-wave 
effects,  while  noticeable,  still  do  not  significantly  alter  the  qualitative  behavior  of  the 
vortex.  Our  choice  to  use  linearized  waves  at  the  amplitudes  chosen  for  the 
primary  cases  is  therefore  a  reasonable  one. 

For  vortices  that  are  linearly  stable  (as  the  one  investigated  here)  and  exhibit  a 
generally  decaying  fluctuation  amplitude,  we  can  expect  another  nonlinear  effect  even 
when  a  is  not  very  large.  For  a  uniform  plane-shear  flow,  the  wrapped-up  filaments 
of  fluctuation  potential  vorticity  are  susceptible  at  long  times  to  secondary  barotropic 
instabilities  associated  with  their  ever-increasing  local  gradient  (Haynes  1987, 
Vanneste  1999).  It  is  likely  that  this  phenomenon  occurs  in  vortices  as  well. 
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MEAN  ENERGY  NORM 


MEAN  ENSTROPHY  NORM 
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Figure  12.  History  of  the  vortex  change  norms  for  (left)  energy,  AEswe(t),  and  (right)  potential  enstrophy, 
AVswe(t),  for  the  primary  cases  with  sustained  external  strain:  CYCLONE  =  solid  line:  ANTICYCLONE  = 
long-dash  line;  QG  =  dash-dot  line;  and  BT  =  dotted  line.  Normalizations  are  in  table  1 . 


However,  its  consequences  should  be  limited  to  the  late-stage  evolution  of  the 
fluctuations  after  the  period  in  which  their  eddy-mean  interaction  is  largely  over, 
so  we  do  not  investigate  that  process  here. 


8.  Sustained  strain 

Section  4  described  how  vortices  react  to  an  impulsive  strain.  Here  we  impose  the  same 
external  strain,  but  after  the  ramp-up  period  ( t  >  tmax/ 2),  the  amplitude  is  fixed 
throughout  the  runs  (i.e.,  F(t)  —  1).  Thus  we  will  see  a  combination  of  both  the 
strain  response  (section  4)  and  the  relaxation  from  that  response  (section  5).  We 
made  runs  for  all  four  primary  cases:  CYCLONE,  ANTICYCLONE,  QG,  and  BT. 

Figure  12  shows  the  mean  vortex  energy  and  enstrophy  change  for  all  the  four  cases. 
Note  the  dramatic  difference  between  CYCLONE  and  ANTICYCLONE:  the  cyclone 
exhibits  a  continual  strong  vortex  weakening,  while  the  anticyclone,  after  an  initial 
weakening  response  (section  4),  begins  to  strengthen  and  continues  to  do  so.  The  QG 
and  BT  vortices  undergo  net  weakening,  with  relatively  small  amplitudes  compared 
to  the  cyclone.  Their  weakening,  however,  occurs  in  an  oscillatory  fashion  that  is  not 
clearly  visible  here  due  to  their  small  amplitude.  In  particular,  the  QG  vortex  has  an 
initial  weakening,  followed  by  a  partial  recovery,  and  then  a  continued  weakening. 
The  BT  vortex  initially  weakens,  but  then  oscillates  until  it  asymptotes  to  a  constant 
negative  value  for  both  mean  energy  and  enstrophy.  This  last  case  is  qualitatively 
very  similar  to  the  results  reported  in  BG99  (their  figure  4). 

As  the  external  strain  turns  on,  all  the  four  cases  initially  have  a  response  as  discussed 
in  section  4.  The  strain  induces  a  wavenumber-two  fluctuation  by  elongating  the 
mean  vortex  into  an  elliptical  shape.  These  fluctuations  then  attempt  to  relax  as 
VRWs  (section  5),  yet  must  do  so  while  still  immersed  in  the  external  strain  flow. 
The  VRWs  begin  to  wind  up  spirally,  but  are  impeded  by  the  fixed  strain  flow. 
This  prevents  axisymmetrization  of  the  vortex,  so  the  vortex  remains  elliptical,  reaching 
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Figure  13.  Swirl  plots  for  qah(x,y)  normalized  by  fV/L  at  1  =  31.4L/K  in  (a)  CYCLONE  (with  contour 
interval  0.117)  and  (b)  ANTICYCLONE  (with  contour  interval  0.025  and  the  sign  reversed)  under  sustained 
external  strain.  Negative  contours  are  dashed,  and  the  zero  contours  are  labeled. 


equilibrium  with  a  standing-wave  fluctuation  pattern  oriented  45°  counter-clockwise 
from  the  initial  strain  induced  location.  This  can  be  seen  in  the  late  time  swirl  plots 
of  normalized  qah  for  both  the  cyclone  and  anticyclone  in  figure  13.  Compare  this 
with  figure  6  which  demonstrates  axisymmetrization.  This  45°  orientation  of  the 
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Figure  14.  Hovmiiller  plots  of  maxx[\(j "  |](r,  t)  under  sustained  external  strain.  The  fields  are  normalized 
by  their  maximum  value  near  r  =  L  and  t  =  3L/V:  (a)  CYCLONE  (max  value  =  0.736  JV/L)  and 
(b)  ANTICYCLONE  (max  value  =  0.095  JVjL).  The  contour  interval  is  0.1. 


vortex  under  sustained  strain  was  also  demonstrated  by  BG99  for  the  BT  case 
(their  figure  7a).  This  is  in  stark  contrast  to  the  behavior  of  an  initially  circular 
vortex  patch  placed  into  a  weak,  steady  straining  flow,  which  oscillates  continually 
between  circular  and  elliptical  states  (Kida  1981;  BG99).  The  reason  for  this  behavior 
is  explained  by  BG99:  a  vortex  patch  does  not  exhibit  the  spiral  wrap-up  of  q  contours, 
thus  arresting  the  irreversible  transfer  of  enstrophy  to  small  scales,  and  making  it 
impossible  for  the  vortex  to  relax  to  the  stable,  steady-state  solution  of  an  elliptical 
vortex  patch  with  its  major  axis  tilted  45°  relative  to  the  strain  axes  (Moore  and 
Saffman  1971). 

As  the  vortex  fluctuations  evolve  into  the  configuration  in  figure  13,  they  interact 
with  the  external  strain  resulting  in  mean  vortex  weakening  for  the  cyclone  and 
strengthening  for  the  anticyclone.  This  can  be  seen  by  identifying  the  dominant  terms 
in  the  vortex  change  azimuthal  velocity  equations  (8)  and  (10)  as 

3*<Vi>  *  -(EAO  »  -(«%■},  (39) 

where  the  fluctuation  relative  vorticity  is 

£  =  Iar(n/)--^M;..  (40) 

However,  the  fluctuations  are  not  lost  through  this  interaction  because  they  are 
continually  reinforced  by  interaction  between  the  mean  vortex  and  the  external  strain 
(even  growing  in  the  outer  region  of  the  anticyclone;  figure  14).  Thus  individual  parcels 
rotate  around  the  vortex,  but  the  waves  are  relatively  locked  in  a  standing  eddy  pattern. 
There  is  some  exception  to  this,  as  seen  by  the  increasing  development  of  fine  structure 
in  the  inner  core  region,  where  there  is  competition  between  the  fixed  external  strain  and 
the  VRW  tendencies  for  spiral  wrap-up  and  radial  propagation. 

These  swirl  plots  highlight  other  interesting  and  telltale  features  of  the  vortex 
evolution  under  sustained  strain.  Both  the  cyclone  and,  to  a  larger  degree,  the 
anticyclone  develop  a  distinct  outer  region,  which  for  the  anticyclone  both  grows 
and  has  the  opposite  sign  to  the  inner  core  region.  The  vortex  strengthening  of  the 
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Figure  15.  Changes  in  (a)  (0)(r)  (normalized  by  m 2  VfL)  and  (b)  (v)(r)  (normalized  by  m2  V)  at  t  =  31.4L/K 
for  CYCLONE  (solid  lines)  and  ANTICYCLONE  (long-dash  lines,  with  the  sign  reversed)  under  sustained 
external  strain. 


anticyclone  (as  compared  with  the  weakening  for  the  cyclone)  is  essentially  due  to 
this  difference.  Looking  back  at  the  impulsive  strain-induced  fluctuation  geopotential 
for  the  different  cases  in  section  4  (figure  4),  the  anticyclonic  vortex  response  in  the 
outer  core  has  a  fluctuation  of  opposite  sign  to  the  inner  core,  while  the  other  cases 
do  not.  Although  the  0(1)  forcing  term  in  (24)  switches  sign  between  cyclones 
and  anticyclones,  the  O(Ro)  terms  (i.e.,  the  second  right-side  group  within  the  square 
brackets  proportional  to  y2)  do  not.  Thus  the  forcing  terms  are  additive  in  cyclones 
and  partially  canceling  in  anticyclones,  such  that  at  large  radius  the  sum  of  the 
anticyclonic  forcing  terms  is  opposite  in  sign  to  the  core.  It  is  this  outer  anomaly 
that  grows  and  propagates  outwards  radially  which  acts  to  strengthen  the  vortex. 
This  can  be  seen  in  the  Hovmuller  plots  of  normalized  \qah\{r,t)  in  figure  14.  These 
demonstrate  not  only  the  radial  propagation  of  the  VRWs,  but  also  their  growth, 
as  well  as  the  development  of  inner  core  fine  structure.  Note  that  the  outer  fluctuations 
do  not  stagnate,  and  that  the  anticyclonic  fluctuation  continues  to  grow  fairly  rapidly  as 
long  as  the  external  strain  is  sustained. 

These  fluctuations  interact  primarily  with  the  external  strain  to  produce  a  mean 
vortex  change  response  in  geopotential  (<j>)  and  azimuthal  velocity  (v).  The  end-time 
values  are  plotted  in  figure  15.  Note  the  correspondence  between  the  outer-region 
location  of  the  fluctuations  and  the  signal  in  (v)  that  in  the  norms  contributes  to 
weakening  for  cyclones  and  strengthening  for  anticyclones.  The  plots  of  (</>)  show  the 
controlling  element  responsible  for  vortex  weakening  in  the  cyclone  and  strengthening 
for  anticyclone  in  the  core  region. 

The  key  difference  between  the  cyclone  and  anticyclone  in  the  finite-Ro  regime  is  the 
existence  of  an  oppositely  signed  potential  vorticity  fluctuation  in  the  outer  core  region 
of  the  anticyclone.  This  allows  for  a  continual  strengthening  of  the  mean  vortex  as 
the  strain-induced  fluctuations  interact  with  the  fixed  external  strain.  This  continual 
strengthening  can  be  seen  even  with  Ro  as  small  as  0.0375  when  Ld/Lw  1,  so  this 
sustained-strain  growth  is  not  just  a  large-Ro  phenomenon.  In  the  cyclone  (and  QG 
and  BT  vortices),  the  outer  and  inner  fluctuations  are  of  the  same  sign.  This  results 
in  the  familiar  vortex  weakening,  demonstrated  earlier  in  this  article,  that  continues 
as  long  as  the  external  strain  is  maintained. 

From  a  meteorological  perspective,  this  anticyclonic  strengthening  mechanism  is 
partly  distinct  from  previous  explanations  of  persistent  blocking  anticyclones  in  the 
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mid-latitude  troposphere.  In  addition  to  viewing  the  blocking  flow  configuration  as  a 
stable,  nonlinear,  stationary  state  (Charney  and  DeVore  1979,  McWilliams  1980, 
Butchart  et  al.  1987),  it  has  been  argued  that  the  anticyclonic  vortex  is  further  strength¬ 
ened  by  up-scale  eddy  vorticity  fluxes  from  remotely  generated  smaller  eddies  that 
impinge  on  the  parent  vortex  (Green  1977,  Shutts  1983).  The  demonstration  here 
that  fluctuations  locally  induced  by  sustained  strain  can  strengthen  an  anticyclone 
with  finite  Ro  suggests  an  alternative  explanation  for  how  atmospheric  blocking  is 
sustained. 


9.  Summary  and  prospects 

Isolated  coherent  vortices  in  the  interior  regions  of  rotating,  stratified  flows  undergo 
cycles  of  disruption  by  other  flows  in  their  neighborhood  and  recovery  by  axisymmetri- 
zation.  In  this  study,  we  have  idealized  this  cycle  to  the  perturbation  of  an  otherwise 
stationary  vortex  in  the  Shallow-water  equations  by  a  transient  external  strain  flow 
followed  by  a  relaxation  phase  in  which  linearized  asymmetric  fluctuations  evolve  as 
Vortex  Rossby  Waves  (VRWs)  by  propagating  in  radius,  wrapping  around  the  mean 
vortex,  decaying  in  energy,  and  forcing  changes  in  the  azimuthally  averaged  vortex. 
The  tendency  of  these  vortex  changes  is  weakening  during  the  straining  phase 
and  strengthening  during  the  recovery  phase,  albeit  typically  only  partially  so  that 
the  outcome  of  the  cycle  is  net  weakening.  Furthermore,  for  intermediate  values  of 
the  vortex  Rossby  and  Froude  numbers,  cyclones  are  more  disruptable  and  have  a 
greater  net  weakening  than  anticyclones  for  reasons  that  are  understandable  in  terms 
of  the  response  to  an  impulsive  strain,  independent  of  the  subsequent  VRW  evolution. 
In  the  case  of  a  sustained  external  strain  flow,  the  same  kind  of  disparity  between 
cyclonic  and  anticyclonic  responses  occurs,  except  that  a  sufficiently  strong  anticyclone 
can  manifest  net  strengthening. 

We  believe  this  is  the  likely  explanation  for  the  many  reported  examples  of  anti¬ 
cyclonic  vortex  dominance  in  the  theoretical  literature  (section  1),  as  well  as  the  appar¬ 
ently  greater  robustness  of  strong  anticyclones  in  nature  away  from  boundaries.  Thus, 
we  believe  that  strain-weakening,  VRW-strengthening,  and  anticyclonic-dominance  are 
generic  behaviors  for  coherent,  interior,  rotating,  stratified  vortices,  although  it  will  be 
necessary  to  explore  well  beyond  the  idealizations  in  this  study  to  adequately  test  this 
hypothesis. 
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